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The dynamics of Brownian motion has widespread applications extending from 

_L | transport in designed micro-channels up to its prominent role for inducing transport 

in molecular motors and Brownian motors. Here, Brownian transport is studied in 
micro-sized, two dimensional periodic channels, exhibiting periodically varying cross 
sections. The particles in addition are subjected to an external force acting along- 
side the direction of the longitudinal channel axis. For a fixed channel geometry, 
the dynamics of the two dimensional problem is characterized by a single dimen- 
sionless parameter which is proportional to the ratio of the applied force and the 
temperature of the particle environment. In such structures entropic effects may 
play a dominant role. Under certain conditions the two dimensional dynamics can 

_J , be approximated by an effective one dimensional motion of the particle in the longi- 

tudinal direction. The Langevin equation describing this reduced, one dimensional 
process is of the type of the Fick-Jacobs equation. It contains an entropic poten- 
tial determined by the varying extension of the eliminated channel direction, and 
a correction to the diffusion constant that introduces a space dependent diffusion. 

/\ • Different forms of this correction term have been suggested before, which we here 

compare for a particular class of models. We analyze the regime of validity of the 
Fick-Jacobs equation, both by means of analytical estimates and the comparisons 
with numerical results for the full two dimensional stochastic dynamics. For the 
nonlinear mobility we find a temperature dependence which is opposite to that 
known for particle transport in periodic potentials. The influence of entropic effects 
is discussed for both, the nonlinear mobility and the effective diffusion constant. 
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1 Introduction 



The phenomenon of entropic transport is ubiquitous in biological cells, ion 
channels, nano-porous materials, zeolites and microfluidic devices etched with 
grooves and chambers. Instead of diffusing freely in the host liquid phase the 
Brownian particles frequently undergo a constrained motion [1,2 3 4 5 6 7 8 9. 10|I11|I12] 
The geometric restrictions to the system's dynamics results in entropic barri- 
ers and regulate the transport of particles yielding important effects exhibiting 
peculiar properties. The results have prominent implications in processes such 
as catalysis, osmosis and particle separation [1.23.45.67.89. 10)1 1 1 II 1 2 j and, as 
well, for the noise-induced transport in periodic potential landscapes that lack 
reflection symmetry (Brownian ratchet systems) [T3|14lll5j or Brownian motor 
transport occurring in arrays of periodically arranged asymmetric obstacles, 
termed "entropic" ratchet devices [IB] . Motion in these systems can be in- 
duced by imposing different concentrations at the ends of the channel, or by 
the presence of external driving forces supplying the particles with the energy 
necessary to proceed. The study of the kinetics of the entropic transport, the 
properties of transport coefficients in far from equilibrium situations and the 
possibility for transport control mechanisms are pertinent objectives in the 
dynamical characterization of those systems. 

Because the role of inertia for the motion of the particles through these struc- 
tures can typically be neglected the Brownian dynamics can safely be analyzed 
by solving the Smoluchowski equation in the domain defined by the available 
free space upon imposing the appropriate boundary conditions. Whereas this 
method has been very successful when the boundaries of the system possess a 
rectangular shape, the challenge to solve the boundary value problem in the 
case of nontrivial, corrugated domains represents a difficult task. A way to cir- 
cumvent this difficulty consists in coarsening the description by reducing the 
dimensionality of the system, keeping only the main direction of transport, 
but taking into account the physically available space by means of an entropic 
potential. The resulting kinetic equation for the probability distribution, the 
so called Fick- Jacobs (FJ) equation, is similar in form to the Smoluchowski 
equation, but now contains an entropic term. The entropic nature of this term 
leads to a genuine dynamics which is distinctly different from that observed 
when the potential is of energetic origin [17] . It has been shown that the 
FJ equation can provide a very accurate description of entropic transport in 
channels of varying cross-section [T7|18|19j . However, the derivation of the 
FJ equation entails a tacit approximation: The particle distribution in the 
transverse direction is assumed to equilibrate much faster than in the main 
(unconstrained) direction of transport. This equilibration justifies the coars- 
ening of the description leading in turn to a simplification of the dynamics, 
but raises the question about its validity when an external force is applied. To 
establish the validity criterion of a FJ description for such biased diffusion in 




Fig. 1. Schematic diagram of a channel confining the motion of forced Brownian 
particles. The half-width uj is a periodic function of x with periodicity L. 

confined media is, due to the ubiquity of this situation, a subject of primary 
importance. 

Our objective with this work is to investigate in greater detail the FJ-approxi- 
mation for biased diffusion and to set up a corresponding criterion describing 
its regime of validity. We will analyze the biased movement of Brownian par- 
ticles in 2D periodic channels of varying cross section and formulate different 
criteria for the validity of such a F J-description. On the basis of our numerical 
and analytical results we recapitulate the striking and sometimes counterintu- 
itive features [17], which arises from entropic transport and which are different 
from those observed in the more familiar case with energetic, metastable land- 
scapes 



2 Diffusion in confined systems 



Transport through pores or channels (like the one depicted in Fig. [T]) may 
be caused by different particle concentrations maintained at the ends of the 
channel, or by the application of external forces acting on the particles. Here 
we will exclusively consider the case of force driven transport. The external 
driving force is denoted by F = Fe x . It points into the direction of the channel 
axis. In general, the dynamics of a suspended Brownian particle is overdamped 
[2T] and well described by the Langevin equation, 



dx -> I — : — — ->,-. 



;i) 



where t is time, x the position vector of the particle, rj its friction coefficient, 
k-Q the Boltzmann constant and T the temperature. The thermal fluctuating 
forces which the surrounding fluid exerts on the particle are modeled by zero- 
mean Gaussian white noise £(£), obeying the fluctuation-dissipation relation 
(&(*) £#)) = 2 Sij 5(1 - i') for z, j = x, y, z. 

In addition to Eq. (CD), the full problem is set up by imposing reflecting bound- 



ary conditions at the channel walls. The form of the channel will be specified 
below. 

To further simplify the treatment of this model we introduce dimensionless 
variables. We measure all lengths in units of the period L, i.e. 

x = Lf, (2) 

where r is the dimensionless position vector of the particle. As unit of time 
t we choose twice the time it takes for the particle to diffusively cover the 
distance L which is given by r = L 2 ^/(fcsT), hence, 

i=rt. (3) 

In these dimensionless variables the Langevin equation reads 

df -> -> 

* = /+«<*)■ (4) 

— * — * — * 

where (£(£)) = and (£(£)£(£')) = 25ij5(t — t') for i,j — x,y,z and where the 
dimensionless force 

f = fe x andf = —— (5) 

k B T 

contains the dimensionless parameter / that characterizes the force as the 

ratio of the work which it performs on the particle along a distance of the 

length of the period and the thermal energy. The corresponding Fokker-Planck 

equation for the time evolution of the probability distribution P(r, t) takes the 

form [221123] : 

^^ = -V-J(f,t), (6a) 

where J(r,t) is the probability current: 

J(r,t) = (f-V)P(r,t), (6b) 



Note that for channels with similar geometry, which are related by a scale 
transformation x — > Xx, A > 0, the transport properties are determined by 
the single dimensionless parameter / which subsumes the period lengths, the 
force and the temperature of the surrounding fluid. 

The reflection of particles at the channel walls leads to a vanishing probability 
current at the boundaries. Therefore, the boundary conditions at the channel 
walls are: 

J(f, t) ■ n = r G channel wall. (7) 

where ft denotes the normal vector field at the channel walls. 



The boundary of a 2D periodic channel which is mirror symmetric about its 
axis is given by the periodic functions y = ±u(x), i.e. u(x + 1) = u>(x) for all 
x, where x and y are the Cartesian components of f. In this case, the boundary 
condition becomes 



do; (a;) 

dx 



fP(x,y,t) — 



dP ^V,t) = (8) 

dy 



at y — ±u>(x). Except for a straight channel with u = const, there are no 
periodic channel shapes for which an exact analytical solution of the Fokker- 
Planck equation ([S]) with boundary conditions (jSj) is known. Approximate 
solutions though can be obtained on the basis of an one dimensional diffusion 
problem in an effective potential. Narrow channel openings, which act as ge- 
ometric hindrances in the full model, show up as entropic barriers in this one 
dimensional approximation [17fl8ll24ll25f2(j||27j . This approach is valid under 
conditions that will be discussed below in some detail. 



3 The Fick-Jacobs approximation 



In the absence of an external force, i.e. for / = 0, it was shown 
that the dynamics of particles in confined structures (such as that of Fig. [1]) can 
be described approximatively by the FJ equation, with a spatially dependent 
diffusion coefficient: 

^-»^ Ml) »f). (9) 



dt dx V dxuj(x) 



This ID equation is obtained from the full 2D Smoluchowski equation upon 
the elimination of the transversal y coordinate assuming fast equilibration in 
the transversal channel direction. Here P{x,t) = J^tLdy P(x,y,t) denotes 
the marginal probability density along the axis of the channel. We note that 
for three dimensional channels an analogue approximate Fokker-Planck equa- 
tion holds in which the function u(x) is to be replaced by ttuj 2 (x) (area of 
cross-section). In the original work by Jacobs [24] the ID diffusion coefficient 
D(x) is constant and equals the bare diffusion constant which is unity in the 
present dimensionless variables. Later, Zwanzig [25] and Reguera and Rubi [26] 
proposed different spatially dependent forms of the ID diffusion coefficient. 



3.1 Spatially dependent ID diffusion coefficients 



The ID diffusion coefficient suggested by Zwanzig results from a systematic 
expansion in terms of the gradient of the boundary function uj{x). In leading 



order he obtained D(x) = 1 — 7u/(x) 2 + . . . , where 7 = 1/3 for the considered 
2D structure (for a 3D structure, 7 = 1/2) and the prime denotes the derivative 
with respect to x. Interpreting this result as the first two terms of a geometric 
series Zwanzig proposed his resumed expression for the ID diffusion coefficient 
reading 

D z (x) = — \ (10) 

1 + 'JUJ'{X) Z 

Reguera and Rubi [26J put forward a different form of the ID diffusion coeffi- 
cient: 

D ™ {x) = (1 + J(s)2yr (11) 

which also can be considered as a re-summation of Zwanzig's perturbational 
result. Yet another form of the ID diffusion coefficient was proposed by Kalinay 
and Percus [27J. For the biased diffusion discussed here the results obtained 
with the Kalinay-Percus diffusion coefficient differ only little from those ob- 
tained by the Reguera- Rubi diffusion coefficient. Therefore we do not further 
consider the Kalinay-Percus diffusion coefficient. 



3.2 Constant bias along the channel direction 



In the presence of a constant force F along the direction of the channel the 
FJ equation (JSJ can be recast into the form 




'121 



with the dimensionless free energy A(x) := E — S = — fx — lnuj(x). In physical 
dimensions the energy is E = ksTE = —Fx (x = xL) and the dimensional 
entropic contribution is S = ksTS = k-^T In a;. For a periodic channel this free 
energy assumes the form of a tilted periodic potential. In the absence of a force 
the free energy is purely entropic and Eq. ffT2|) reduces to the FJ equation (jUJ). 
On the other hand, for a straight channel the entropic contribution vanishes 
and the particle is solely driven by the external force. 



3.3 Nonlinear mobility and effective diffusion 



Key quantities of particle transport through periodic channels are the average 
particle current, or equivalently the nonlinear mobility, and the effective dif- 



fusion coefficient. For a particle moving in a one dimensional tilted periodic 
potential the heights AV^ of the barriers separating the potential wells provide 
an additional energy scale apart from the work of the force FL and the thermal 
energy ksT. Hence, at least two dimensionless parameters, say AE / '(fe^T) and 
FL/(kBT) govern the transport properties of these systems. In contrast, as 
already noted in the context of the full 2D model the transport through chan- 
nels is governed by the single dimensionless parameter / = FL/(ksT). This, 
of course, remains to hold true in the one dimensional approximation which 
models the transversal spatial variation in terms of an entropic potential. 

For any non negative force the average particle current in periodic structures 
can be obtained from 1281: 



(x) = (t{x ^x + I))' 1 (13) 

where t(a —>■ b) denotes the first time of a particle which starts at x = a to 
arrive at x — b. The angular brackets refer to an average over the fluctuating 
force. Note that the resulting mean first passage time (t(xo — > x + l)) diverges 
for a vanishing force and consequently leads to a vanishing current. A positive 
force prevents the particle to make far excursions to the left, hence leading 
to a finite mean first passage time as well as a finite current. Within the one 
dimensional approximation, cf. Eq. (1121) . the moments of the first passage time 
t(a — > b) can be determined recursively by means of 

(t n ( a -►&))= n / dx —— exp (A(x)) f dy exp (-A(y)) 

Ja D[X) J-oo 

(t n -\y^b)). (14) 
For n = the starting value of the iteration is given by (t°(a — > b)) = 1. 
The nonlinear mobility fi(f) is defined by 

M/) = y- ( 15 ) 

Using eqs. (TT3|) and (fi~4l one can obtain the following Stratonovich formula 
for the nonlinear mobility [17J 

l-exp(-J) 

M/) = — 71 > ( 16a ) 



/ fdzl{zj) 
Jo 



where 

h-\z) 



Viz) Jz-l 



D cS == JU : T -3 , (19a) 



depends on the dimensionless position z, the force / and the shape of the tube 
given in terms of the half width u(x) and its first derivative. 

The effective diffusion coefficient is defined as the asymptotic behavior of the 
variance of the position 

n , (s 2 (*)) ~ (*(*)) 2 nn 

^ cff = t hm . (17) 

It is related to the first two moments of the first passage time by the expression 
[28JT29] : 

n (t 2 (xo^xo + l))-(t(x ^xo + l)} 2 

UeS ~ ^777 r~m3 • ^° 

2 (t(x -> x + i)y 

After some algebra it can be transformed to read 

f\z f dz/V(z,z,f) 

JO Jz-1 

f'dzlizj) 

Jo 

where 

JVC*. *. /) : = fgf ^y ['& /)] 2 exp(-/ z + fz). (19b) 

5.^ Exact numerics for the 2D channel 



The predicted dependence of the average particle current and the effective 
diffusion coefficient, predicted above, was compared with Brownian dynamic 
simulations performed by a numerical integration of the Langevin equation 
Eq. (jl]), within the stochastic Euler-algorithm. The shape of the exemplarily 
taken 2D channel is described by 

lu(x) := a sin(27ra;) + b , (20) 

where b > a. The sum and difference of the two parameters a + b and b — a 
give half of the maximal and the minimal width of the channel, respectively. 
Moreover, a controls the slope of the channel walls which determines the one 
dimensional diffusion coefficient D(x). 

For the considered channel configuration, cf. Eq. (|20|) . the boundary condition 
becomes u(x) = a (sin(27rx) + k), where k = b/a = 1.02 throughout this 
paper. For a we chose values between 1 and l/2ir. This choice of parameters 
corresponds to rather short channels for a = 1 and a more elongated one for 
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Fig. 2. (color online) The dependence of the nonlinear mobility of a Brownian par- 
ticle within a 2D channel with half width uj(x) = sin(2-7T2;) + 1.02 is depicted versus 
the dimensionless force / = FL/k^T. The symbols correspond to the numerically 
simulated exact values, cf. Eq. (|2ip . and the lines to the analytically calculated, ap- 
proximative values, cf. Eq. ([TBI) , with the diffusion coefficients: Dq = 1 (solid line); 
Drr(x) (dashed line), cf. Eq. (fTTj) . and Dz(x) (dashed-dotted line), cf. Eq. ([TO]) . 



a = 1/27T. In all cases the width of the widest opening of the channel is larger 
by a factor of 100 than the width at narrowest opening. One may therefore 
expect strong entropic effects for these channels. The particle current and 
effective diffusion coefficient were derived from an ensemble- average of about 
3 • 10 4 trajectories: 



(x) 



lim 

t— >oo 



<*(*)> 



(21) 



and Eq. i fPTj) . respectively. 



Fig. [2] demonstrates that the resumed one dimensional diffusion coefficient 
D-rr[x) leads to a considerably better agreement with the numerical results 
than the constant diffusion and the diffusion coefficient Dz(x) proposed by 
Zwanzig. Therefore we used Dn R (x) for all following calculations. 



Validity of the Fick-Jacobs description in the presence of a con- 
stant bias 



The reduction of dimensionality leading to the FJ equation relies on the as- 
sumption of equilibration in the transverse direction which results in an al- 
most uniform distribution of the transversal positions y at fixed values of the 
longitudinal coordinate x. One can formulate two different sets of criteria de- 
termining first, whether the FJ equation describes the relaxation towards the 
stationary state, say in presence of periodic boundary conditions, or second, 
whether the stationary state, but not necessarily the relaxation towards this 
state, can be described by the FJ equation. 



Although we here are mainly interested in the second, weaker, type of criteria 
which are sufficient to guarantee the validity of the transport properties pre- 
dicted by the FJ equation, we shortly formulate stronger criteria which must 
be satisfied if the F J equation is employed to determine the relaxation towards 
equilibrium. Then, of course, the time scale of equilibration in transversal di- 
rection tt must be short compared to the relevant time scales in longitudinal 
direction, tl [18]. The time scale tt can be estimated by the time to diffusively 
cover the widest transversal distance of the channel. It therefore is given by 

r T = 2(a + bf. (22) 

The time scales characterizing the longitudinal motion are the diffusion time 
t al over the length of a period, which is one in our dimensionless variables, 
and the time t {L it takes to drag a particle over this distance by applying the 
force /. These times are given by 

r dL = 1/2 and r iL = 1/f (23) 

Hence, a necessary condition, that the FJ equation reliably describes transient 
processes can be formulated as 

2(a + 6) 2 <min(l/2,l//) (24) 

This condition is fulfilled only for rather elongated channels being at least five 
times as long as wide. 

As already mentioned in the presence of periodic boundary conditions in the 
longitudinal direction, the particle distribution described by the Fokker-Planck 
equation (jBj) approaches a stationary distribution. Even in the case if the first 
criterion (124"]) is violated the stationary solution of the FJ equation (Tl2|) may 
still yield the correct marginal probability density provided transversal cuts 
of the two dimensional stationary probability density are practically constant. 
Such a uniform distribution in transversal direction strictly holds in the ab- 
sence of externally imposed concentration differences if the force / vanishes 
or if the channel is straight. For channels with varying width the narrow po- 
sitions confine the positions of particles. From there they are dragged by the 
force and at the same time they perform a diffusive motion until the channel 
narrows again. The required uniform distribution in the transversal direction 
can only be achieved if the diffusional motion is fast enough in comparison to 
the deterministic drift under the influence of /. In other words, the diffusional 
spreading within the time the force drags the particle from the narrowest to 
the widest place in the channel must be at least of the order of the widest 
channel width. This leads to the second, weaker, criterion 

< a+6 > 2 ^ * f^hwrw*^ (25) 
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Fig. 3. The dependence of the critical value of the dimensionless force f c on 
the parameter a for 2D channels defined by the dimensionless boundary function 
uj{x) = a (sin(27rx) + 1.02) is depicted. For f < f c the accuracy of the Fick-Jacobs 
approximation is ~ 1% (in comparison with the exact 2D simulational results). 
The solid line demonstrates the a _2 -dependence of the critical values predicted by 
Eq. (|25p . The inset depicts the same data on a logarithmic scale. 



Eq. (125]) provides an estimate of the minimum forcing above which the FJ 
description is expected to fail in providing an accurate description of the 
transport properties in the long time limit. The quantitative value of the 
critical force depends on the level of the prescribed accuracy. The criterion 
demonstrates how the validity of the equilibrium approximation depends on 
the relevant parameters of the problem and is concordant with that found for 
a different scaling in Ref. [T8] . 



In order to test the accuracy of the FJ description, we evaluated the behavior 
of the nonlinear mobility as a function of the scaled force /, for different values 
of a according to Eq. ( TTBT) and compared it with numerical simulations of the 
corresponding full two dimensional problem. The value of the dimensionless 
force / up to which the FJ approximation with spatially dependent diffusion 
coefficient D RR (x) provides an accurate description depends on a just as pre- 
dicted by Eq. (1251) . cf. Fig. [3j For large values of a the FJ equation starts to 
deviate from the numerically exact behavior already for rather small forces /, 
whereas for small values of a larger forces may be applied without violating 
the FJ equation. 

In Fig. H] the dependence of the nonlinear mobility on the force / is displayed 
for two different values of the parameter a characterizing the geometry of the 
channel. For the channel with a = 1 which, at its widest opening, is approx- 
imately four time as wide as it is long, the predictions of the equilibration 
assumptions fails at smaller /-values than in the case of a = l/2ir. 
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Fig. 4. (color online) The numerically simulated (symbols) and analytically calcu- 
lated (cf. Eq. (fl~6j) - lines) dependence of the scaled nonlinear mobility /u(/) vs. the 
dimensionless force / = FL/k-^T is depicted for two 2D channel geometries. For 
both channels the scaled half-width is given by to{x) = a (sin(27rx) + 1.02); a = 1 : 
diamonds and solid line(blue), a = 1/(2tt) : circles and dashed line (green). The 
dotted line indicates the deterministic limit //(/) = (x)/f = 1. 

5 Transport characteristics: Anomalous temperature dependence 
and enhancement of diffusion 



Transport in one dimensional periodic energetic potentials behaves very dif- 
ferently from one dimensional periodic systems with entropic barriers. The 
fundamental difference lies in the temperature dependence of these models. 
Decreasing temperature in an energetic periodic potential decreases the tran- 
sition rates from one period to the neighboring by decreasing the Arrhenius 
factor exp{— AV/^ksT)} where AV denotes the activation energy necessary to 
proceed by a period [20] • Hence decreasing temperature leads to a decreasing 
mobility. For a one-dimensional periodic system with an entropic potential, a 
decrease of temperature leads to an increase of the dimensionless force param- 
eter / and consequently to an increase of the mobility, cf. Fig. HI 

On the other hand, the dependence of the dynamics on the geometry param- 
eter a clearly reflects the entropic effects on the mobility. A channel with a 
larger a value has wider openings and therefore provides more space where 
the particle can sojourn. This longer residence time within a period of the 
channel diminishes the throughput and consequently the mobility. This is cor- 
roborated by the results of our calculations depicted in Fig. HJ For all values 
of /, an increase in value of a leads to a decrease in the mobility. This holds 
not only where the FJ equation applies but also for large values of / where it 
fails. 



Another interesting effect can be observed for the effective diffusion if looked 
as a function of the force /. Already the expression for the effective diffusion 
(1191) which follows rigorously from the FJ equation displays a maximum as a 
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Fig. 5. (color online) The numerically simulated (symbols) and analytically cal- 
culated (cf. Eq. (fT9|) - lines) dependence of the effective diffusion coefficient D c r 
is depicted vs. the dimensionless force / = FL/k^T for two channels in 2D. For 
both channels the scaled half-width is given by to{x) = a (sin(27rx) + 1.02); a = 1 : 
diamonds and solid line(blue), a = 1/(2-71") : circles and dashed line (green). 

function of / which may even exceed the value 1 of the bare diffusion, cf. Fig. EJ 
For / — > oo the periodic stationary distribution approaches a delta function 
along the x axis and the effective diffusion approaches the bare value 1. If one 
decreases the force to finite but still large values then the stationary distri- 
bution acquires a finite width in the transversal direction with a "crowded" 
region in front of the narrowest place of the channel Ref. [18] . The transport 
becomes more noisy and consequently the effective diffusion exceeds the bare 
value 1. On the other hand if one starts at / = the entropic barriers dimin- 
ish the diffusion such that the effective diffusion is less than bare diffusion. 
Consequently, somewhere in between there must be a value of / with maxi- 
mal effective diffusion. For a = 1 and b = 1.02 the value of the force at the 
maximal effective diffusion is outside the regime of validity of the F J equation. 
The numerical simulations give a much more pronounced peak of the effective 
diffusion. For the less entropic channel with a = 1/27T the maximum is at the 
border of the regime of validity of the F J equation, but the enhancement of the 
effective diffusion constant is less pronounced than for the larger value a = 1. 
These observations lead us to the conclusion that entropic effects increase the 
randomness of transport through a channel and in this way decrease the mo- 
bility and increase the effective diffusion. A similar enhancement of effective 
diffusion was found in titled periodic energetic potentials [28j2 9ll30] . 



6 Conclusions 



In summary, we demonstrated that transport phenomena in periodic channels 
with varying width exhibit some features that are radically different from con- 
ventional transport occurring in energetic periodic potential landscapes. The 
most striking difference between these two physical situations lies in the fact 
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that for a fixed channel geometry the dynamics is completely characterized by 
a single parameter / = FL/(ksT) which combines the external force F caus- 
ing a drift, the period length L of the channel, and the thermal energy ksT, 
which is a measure of the strength of the acting fluctuating forces. Transport 
in periodic energetic potentials depends, at least, on one further parameter 
which is the height of the highest barrier separating neighboring periods. This 
leads to an opposite temperature dependence of the mobility. While the mobil- 
ity of a particle in an energetic potential increases with increasing temperature 
the mobility of a particle in a channel of periodically varying width decreases. 
The incorporation of the spatial variation of the channel width as an entropic 
potential in the FJ equation allows a qualitative understanding of the depen- 
dence of the transport properties on the channel geometry. 

The effective diffusion exhibits a non monotonic dependence versus the di- 
mensionless force /. It starts out at small / with a value that is less than 
the bare diffusion constant, reaches a maximum with increasing / and finally 
approaches the value of the bare diffusion from above. 

It is known from the literature that, under certain conditions, the two-dimensi- 
onal Fokker-Planck equation governing the time dependence of the probability 
density of a particle in the channel can be approximated by a one-dimensional 
Fokker-Planck equation: It is termed the Fick- Jacobs equation and contains 
an entropic potential and a position dependent diffusion coefficient. Various 
forms of the diffusion coefficient can be found in the literature. A comparison 
of different forms for a particular channel geometry leads to the conclusion 
that the expression recently suggested by Reguera and Rubi yields the most 
favorable agreement. In principle the FJ equation describes both the transient 
behavior of a particle and also the stationary behavior of the particle dynamics 
which is approached in the limit of large times, provided appropriate boundary 
conditions confining the motion in the direction of the channel axis are applied. 
In order to study stationary transport, periodic boundary conditions must be 
invoked. We formulated criteria for the validity of the FJ equation for both the 
transient and the stationary regimes and found that the restrictions imposed 
by the criterion in the stationary regime are much less serious than those for 
the transient dynamics. The estimates, which are based on simple dynamical 
arguments, were corroborated by our numerical simulations. 

We restricted our analysis to two dimensional channels. A generalization of 
the presented methods to three dimensional pores with varying cross section 
is in principle straight-forward. We also confined ourselves to channels with a 
mirror symmetry about a vertical axis which in the present case can be chosen 
at x = 1/4. For periodic channel shapes without this symmetry ratchet like 
transport can be expected even if the unbiased force / of vanishing temporal 
average changes periodically in time. 
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